home *** CD-ROM | disk | FTP | other *** search
/ IRIX 6.2 Development Libraries / SGI IRIX 6.2 Development Libraries.iso / dist / complib.idb / usr / share / catman / p_man / cat3 / complib / dtpcon.z / dtpcon
Text File  |  1996-03-14  |  3KB  |  133 lines
  1.  
  2.  
  3.  
  4. DDDDTTTTPPPPCCCCOOOONNNN((((3333FFFF))))                                                          DDDDTTTTPPPPCCCCOOOONNNN((((3333FFFF))))
  5.  
  6.  
  7.  
  8. NNNNAAAAMMMMEEEE
  9.      DTPCON - estimate the reciprocal of the condition number of a packed
  10.      triangular matrix A, in either the 1-norm or the infinity-norm
  11.  
  12. SSSSYYYYNNNNOOOOPPPPSSSSIIIISSSS
  13.      SUBROUTINE DTPCON( NORM, UPLO, DIAG, N, AP, RCOND, WORK, IWORK, INFO )
  14.  
  15.          CHARACTER      DIAG, NORM, UPLO
  16.  
  17.          INTEGER        INFO, N
  18.  
  19.          DOUBLE         PRECISION RCOND
  20.  
  21.          INTEGER        IWORK( * )
  22.  
  23.          DOUBLE         PRECISION AP( * ), WORK( * )
  24.  
  25. PPPPUUUURRRRPPPPOOOOSSSSEEEE
  26.      DTPCON estimates the reciprocal of the condition number of a packed
  27.      triangular matrix A, in either the 1-norm or the infinity-norm.
  28.  
  29.      The norm of A is computed and an estimate is obtained for norm(inv(A)),
  30.      then the reciprocal of the condition number is computed as
  31.         RCOND = 1 / ( norm(A) * norm(inv(A)) ).
  32.  
  33.  
  34. AAAARRRRGGGGUUUUMMMMEEEENNNNTTTTSSSS
  35.      NORM    (input) CHARACTER*1
  36.              Specifies whether the 1-norm condition number or the infinity-
  37.              norm condition number is required:
  38.              = '1' or 'O':  1-norm;
  39.              = 'I':         Infinity-norm.
  40.  
  41.      UPLO    (input) CHARACTER*1
  42.              = 'U':  A is upper triangular;
  43.              = 'L':  A is lower triangular.
  44.  
  45.      DIAG    (input) CHARACTER*1
  46.              = 'N':  A is non-unit triangular;
  47.              = 'U':  A is unit triangular.
  48.  
  49.      N       (input) INTEGER
  50.              The order of the matrix A.  N >= 0.
  51.  
  52.      AP      (input) DOUBLE PRECISION array, dimension (N*(N+1)/2)
  53.              The upper or lower triangular matrix A, packed columnwise in a
  54.              linear array.  The j-th column of A is stored in the array AP as
  55.              follows:  if UPLO = 'U', AP(i + (j-1)*j/2) = A(i,j) for 1<=i<=j;
  56.              if UPLO = 'L', AP(i + (j-1)*(2n-j)/2) = A(i,j) for j<=i<=n.  If
  57.              DIAG = 'U', the diagonal elements of A are not referenced and are
  58.              assumed to be 1.
  59.  
  60.  
  61.  
  62.  
  63.                                                                         PPPPaaaaggggeeee 1111
  64.  
  65.  
  66.  
  67.  
  68.  
  69.  
  70. DDDDTTTTPPPPCCCCOOOONNNN((((3333FFFF))))                                                          DDDDTTTTPPPPCCCCOOOONNNN((((3333FFFF))))
  71.  
  72.  
  73.  
  74.      RCOND   (output) DOUBLE PRECISION
  75.              The reciprocal of the condition number of the matrix A, computed
  76.              as RCOND = 1/(norm(A) * norm(inv(A))).
  77.  
  78.      WORK    (workspace) DOUBLE PRECISION array, dimension (3*N)
  79.  
  80.      IWORK   (workspace) INTEGER array, dimension (N)
  81.  
  82.      INFO    (output) INTEGER
  83.              = 0:  successful exit
  84.              < 0:  if INFO = -i, the i-th argument had an illegal value
  85.  
  86.  
  87.  
  88.  
  89.  
  90.  
  91.  
  92.  
  93.  
  94.  
  95.  
  96.  
  97.  
  98.  
  99.  
  100.  
  101.  
  102.  
  103.  
  104.  
  105.  
  106.  
  107.  
  108.  
  109.  
  110.  
  111.  
  112.  
  113.  
  114.  
  115.  
  116.  
  117.  
  118.  
  119.  
  120.  
  121.  
  122.  
  123.  
  124.  
  125.  
  126.  
  127.  
  128.  
  129.                                                                         PPPPaaaaggggeeee 2222
  130.  
  131.  
  132.  
  133.